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Abstract. 

We study the full susceptibility of the Ising model modulo powers of primes. 
We find exact functional equations for the full susceptibility modulo these primes. 
Revisiting some lesser-known results on discrete finite automata, we show that 
these results can be seen as a consequence of the fact that, modulo 2^, one cannot 
distinguish the full susceptibility from some simple diagonals of rational functions 
which reduce to algebraic functions modulo 2^, and, consequently, satisfy exact 
functional equations modulo 2^. We sketch a possible physical interpretation 
of these functional equations modulo 2^ as reductions of a master functional 
equation corresponding to infinite order symmetries such as the isogenies of elliptic 
curves. One relevant example is the Landen transformation which can be seen 
as an exact generator of the Ising model renormalization group. We underline 
the importance of studying a new class of functions corresponding to ratios of 
diagonals of rational functions: they reduce to algebraic functions modulo powers 
of primes and they may have solutions with natural boundaries. 
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1. Introduction 

Despite the enormous progress made over the last 75 years in the study of 
(Yang-Baxter) integrable models in lattice statistical mechanics and enumerative 
combinatorics, there still remain many important unanswered questions. 

One of the most intriguing is the susceptibility of the two-dimensional Ising model. 
The closed form expressioiQ of the partition function was obtained by L. Onsager [2] 


f It can be rewritten in a simpler 4 F 3 hypergeometric form, see [I]. 
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in 1944, and the spontaneous magnetization was obtained a few years later by both 
Onsager (unpublished) and Yang [3]. However, after more than 70 years, a closed 
form expression for the full susceptibility still eludes us. Accordingly, understanding 
the nature of this function remains a challenging problem. 

Forty years ago, Wu, Barouch, McCoy and Tracy [4] showed that the full 
susceptibility of the square-lattice Ising model can be decomposed as the infinite 
sum of holonomic n-fold integrals [3 El H 0 m, denoted In the last decade 

the linear differential operators corresponding to the first up to n = 6, were 

obtained, underlying the role of the elliptic curve parametrization cni, but showing 
also the emergence of (at least) one Calabi-Yau ODE, and beyond, of linear differential 
operators with selected differential Galois groups [IIllTailS]. A complete description 
of the singular points of the linear differential operators corresponding to the first 
few y("^’s has also been obtained 0 nn HU ng. Despite being an infinite sum of 
holonomic n-fold integrals, the full susceptibility is not a holonomic function m- 

Further, in a recent paper it has been shown that these n-fold integrals are 
actually diagonals of rational functions [HUS]. Consequently their series expansions 
are such that modulo any prime, or power of a prime, they can be identified with the 
series expansions of an algebraic function [SI [S] • These properties were explicitly 
shown in the case of in section (3.1) of [S]- In particular it was shown that 

H{w) = (defined in Section 0 see dSJ), satisfies, modulo 2, the quadratic 

equation 

H{wY +w- H{w) -I-= 0 mod 2. (1) 


and modulo 3, the polynomial equation of degree nine 

Pq ■ H{w)^ + ■ ps- H{w)^ + ■ Pi ■ H{w) + = 0, (2) 


where: 


Pq^^ = + w'^ — w'^ — w + 1, 

p(2) = ^37 _ ^36 ^ ^35 _ ^33 ^ ^31 _ ^30 


,,28 




W 


..22 


_ ^21 _ ^18 _ ^16 


\13 


,^14 _^12 _^10 


Pi = (I - w) 

pg = (li; + 1)3 (^2 ^ 


w — w° — — I, 

15 ,'.,.4 ...2 


P3 = (i„4 ^ 1)48 (1 _ u;)n> (^4 _ 


\24 


( 3 ) 


Since all the x^”^’s are diagonals of rational functions [Sj: similar results are expected 
for any x^"^ modulo any prime p, and, beyond, modulo any power of a prime p'' . As 
a consequence of the Fermat relations, = a, modulo p, one can expect relations, 
like ([U or ([2|), to be expressible as functional equations where H[wY is replaced by 
H{w^). Now, the full susceptibility x is not the diagonal of a rational function, indeed 
it is not even holonomic [suni- Therefore, for the full susceptibility, one cannot 
expect relations like © or (0) to exist. Due to the complexity of this functiorjH, one 
might not expect, at first sight, such functional equations for the full susceptibility. 

However, as we show below, the full susceptibility, when expressed in the 
appropriate expansion variable, does satisfy some surprisingly simple functional 
equations modulo certain primes, or power of primes. 

These exact results show that the full susceptibility reduces to an algebraic 
function, modulo certain primes, or powers of primes, and thus sheds new light 
on the integrable character of this very important function in physics. We consider 


I Which has, for instance, a natural boundary 115111611771 . 
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this a surprising result: we certainly did not expect such simple results for the full 
susceptibility. This gives us considerable incentive to systematically study other non- 
holonomic physical series modulo primes or powers of primes. It will be interesting 
to see whether this is an exceptional result, in which case it sheds more light on the 
susceptibility, or a common occurrence, in which case we need to explain why. 

2. Definitions and some known results on the full susceptibility. 

In 1976, Wu, McCoy, Tracy and Barouch [4] showed that the susceptibility could be 
expressed as an infinite sum of contributions, known as n-particle contributions 
The low-temperature series were given by the case n even, and the high-temperature 
series by n odd. More precisely the low temperature susceptibility is given by m 

kT-xdw) = -jcdw) = (1 - l/s'^)^ • 

= ( 4 ) 

in terms of the self-dual temperature variable w = ^s/(l-|-s^), where s = 

sinh(2 J/A:T), and sl = 1/sinh(2J/fcT). The high temperature susceptibility is given 
by m 

kT-xniw) = - ■ (1-s^)i ■ xniw) = ^ • (1 - • V (5) 

s s ^^ 

Remarkably loM series expansions with respectively 2042 and 2043 coefficient^, 
have been obtaineqS [20] for xl{w) and xh{w), namely 

Xl{w) = 4i(;'‘ -f SOw® -b 1400 w® -b 23520-b 388080 -b 6342336 

-b 103062976 -b 1668639424 -b 26948549680 + ■ ■ ■ + c[oIq w^°®® -b • • • (6) 

and 

XHiw) = 2w-b 8 -b 32 w® -b 128 -b 512 w® -b 2048 w® -b 8192 -b 32768 w® 

-b 131080 -b 524288 + 2097440 -b • • • -b 4o43 + ■■■ (7) 

It is worth comparing these two series with the series corresponding to the first 
x(")(w) in the two infinite sums (gj) and ([S]), namely : 

Xl\w) = 4^“^ • 2i^i([^, ^],[3], 16 u>2) 

= -b80w® -bl400w® -b 23520 -b 388080 -b 6342336 

-b 103062960 w^® -b 1668638400 w®® -b 26948510160 + ■■■ (8) 

and 

yWiw) = ^ = 2w +8w^ +32w^ +128w^ -b512w® -b 2048 w® -b 8192 

-b 32768 w® -b 131072 -b 524288 -b 2097152 -b • • • 

It is known that = 0{w^^), so that the coefficients are the same up to for 
the low-temperature series, and up to w® for the high-temperature series. Further, 

f The low temperature series XZ/(^)> being an even function, means that the expansion is known up 
to the coefficient of 

:( Using an algorithm adapted from the Fortran algorithm in El- 
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one observes that the ratio of the coefficients for xl and x^^ (resp. xh and is 
very close to 1. 

The series expansion for (w) readsQ 


^ 64 + 2470 w^° + 74724 + 1954688 + 46428552 

+ 1029903288^2® + 21716367896 u;®° + 440440693418 + 8663350828976 

+ 166258457615526 w®® + 3126949985578700 w®® + 57833406662680980 
+ 1054656431047823680 ^"*2 + 19003412267837223432 + • • • (10) 


The series expansion for Xl^ (w) reads 
~( 6 )/ '> 

= re®® +144 w®® + 11306 u>^° +641604^^2 .^ 29455804 

26 

+ 1161654484 w^® + 40827303872 w'*® + 1310513628660 w®° 

+ 39090651539936 w®2 + 1097452668063296 w®"* + 29281457807054052 w®® 

+ 748130523334531340 w®® + 18414177309344582452 w®° + ••• (11) 

The difference between xl{w) and Xl \^) reads: 

- y^2) ^ 16 w^® + 1024 w^® +39520 ri;20 + 1195 5 84 y;22 3^275008^24 

+ 742856832 ri;26 + 16478452608 w;28 + 347461886336 m;®° + 7047051094688 w;® 2 
+ 138613613263616 w®'^ + 2660135321848480 w®® + ••• (12) 

The difference between xl{w) and X^^\w) + X^^\w) reads: 

XL - (Xl'’+ Xl^) = 64w®® +9216w®® + 723584ri;40 + 4 1 06 2 6 56^"^2 

+ 1885171456 +74345886976 w^®+2612947447808 w^® + ••• (13) 

The difference between xl{w) and X^^\w) + X^^\w) +x^^\w) reads: 

XL - {Xl'’ + Xl'’ + X^l ) = 256w®^ + 65536w®® + 8815104w®® 

+ 829038592 w™ + 61219149824 ri;^2 + 3779726083072 w^4 
+ 202925982372864 w^® + 9729999547422720 w^® + 424756293921653248 w®® 
+ 17127494149322319872 w;®^ + 645117850681779326976 w®^ + • • • (14) 


The series expansion for (w) reads 

^ = 1^9 + 36 lyii + 4 + 384 ^13 j^gg ^14 j^3532 ^i5 gQ84 ^le 

+ 357391 + 153484 w^® + 6556516 + 3440964 + 116449960 

+ 71553656^22 + 2022814844 u>23 + 1413292572 + 34583048616^25 

+ 26900157072 ?«26 + 584324509812 + 498048104276 u>2® + • • • (15) 


f Since there is an overall integer of the form 2"' for all the coefficients of the or X^^ (w) 

series, we divide them, in the following, by an appropriate power of 2^ factor. The series expansion 
remains an expansion with (smaller) integer coefficients. 
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The series expansion for (w) reads 
~(5)^ N 

+ 100 + 5652 + 4 + 238032 + 484 + 8323743 

+ 32436 + 255716632 ^35 + 1592488^36 + 7139250236 + 63994900 

+ 185181953320^39 + 2231760988+ 4531508893397 + 69986224204 

+ 105775797597812 ^43 + 2020409460692^34 +2374723605151320^45 
+ 54584651129624^46 + 51602310149637388 m; 47 + 1396760803374712 ^48 
+ 1090696414153653447 u;49 + ••• (16) 

Comparing xh(w) and the sum (w) + X^^('U') one finds that these two series 
are the same up to 0{w'^^), as expected: 

XH - (Xff + Xff ) = 32^95 + 3200^97 + 180864 ^99 +128^3° + 7617024 u;34 

+ 15488^32 +266359776 u;33 + 103 7 9 52 ri;34 +8182932224 u;35 + ... (17) 


XH - (Xff + Xff + Xff ) = 128 u;49 + 25088 + 2621952 ^63 

+ 194185216 u;65 +512^66 + 11431676800^67 + 115200^68 

+ 569065324032 ri;69 + 13 7 0 9 8 24 ri;60 + ... ( 18 ) 


Since the modulus of elliptic functions parametrising the Ising model [iniiis] is 
k = s^, with the conjectured natural boundary [16] corresponding to the unit k or 
s circle, it is natural to introduce series expansions in the s or sl variables. In 
fact, we have studied series expansions in the v = sl/2 = 1/(2 s) variable in the 
low-temperature regime, and the v = s/2 variable in the high-temperature regime, 
in order to have series with integer coefficients (instead of rational coefficients with 
denominators of the form 2"). The corresponding low and high temperatures series 
Xl{v) = xl{sl/2.) and xh{v) = xh{s/2) read respectiveljlil 
Xl{v) = 4?;4 +16^6 +104n® +416^4° + 2224?;42 + 8896 ^ 44 + 43 8 40 ^46 

+ 175296 ?;48 + 82 5 6 48 + 3300480 + ... + cf^^^ + • • • (19) 

and 


Xh(v) = 1 +4?; +12z;9 +32v^ +76v'^ +176^6 +400^6 + 896 + 1964 u® 

+ 4256?;® +9184?;40 + 197 28 ?;44 + 41952 ^4® + ••• + + ••• (20) 

which can be compared with Xl^ = (1 “ 16?;4)4/4 . 




43 

= 4?;4 


16 n® + 104?;* 
+ 175296 ?;48 


.1 +4?;2/ "V*2’ 2*’* 4’ Vl +4?;2. 

+ 416?;40 + 2224?;42 + 8896 ?;44 + 43 8 24 ?;46 
825104 ?;®° +3300416 ?; 2 ® + ••• ( 21 ) 


and x^H = (1 - s'‘)4/4/s • x^^ 


Xff 


(1 — s4jl/4 g / 1 — 16?;4 \ 1/4 

S (1 — s)2 +1 — 2v)^ ) 

+ iv +12?;® +32?;3 +76?;4 +176?;6 +400?;° + 896 ?;’'+ 1960 ?;® 
+ 4256 ?;® + 9184 ?;4° + 19712 ?;44 + 41888 ?;4® + • • • (22) 


t Throughout this paper the x functions of the variable w, while the x functions of the 
variable v. 
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As must be the case, the coefficients are the same up to for the low-temperature 
series, and up to u® for the high-temperature series, and, beyond, the ratio of the 
coefficients for xl(v) and (resp. xh(v) and x^h very close to 1. For 

/o'V 

the low-temperature series expansion the difference between Xl{p) and Xl \^) reads: 

Xl{v) - = 16^1® -h 544^2° +64v^^ -H 13056-^2944^^® -h 272512 

-h 88448 -h 5286560 -h 2201856 -H 98136096 u®® -h • • • (23) 

It is worth recalling that the very long (low and high temperature) series 
expansions have been obtained for the full susceptibility as a consequence of 
a quadratic finite difference Painleve functional equation yielding an 

polynomial algorithm for calculating the series. This series is therefore ^^algorithmically 
integrable’'. Furthermore the n-fold integrals of the infinite sum decomposition [1] , the 
x(")’ 3 , have been shown to be highly selected holonomic functions, namely diagonals 
of rational functions m- 

These properties (“algorithmic integrability”, infinite sums of diagonals of rational 
functions, ...) suggest that transcendental non-holonomic functions such as the full 
susceptibility of the square Ising model, should correspond to a “rather special class” 
of non-holonomic functions, which require new concepts and tools to characterize and 
analyze them. 

Obtaining such remarkably long series for the full susceptibility was a 
computational “tour de force,” and it is likely that these series have much more to 
tell us. To date they have only been used to obtain some results on x*'^^ and x*'^^ 
to confirm exact results [Til ITS] on the singularities of the linear ODEs of the x*'"^’s, 
and to clarify the natural boundary scenario |16j . 

In the following sections we revisit these remarkably long series from a new finite 
automaton |21] perspective, which in effect means considering the various series modulo 
various integers, in particular, taking a “p-adic” perspective [22], modulo integers that 
are integer powers of primes. 

3. Functional equations modulo 2” for the full susceptibility. 

3.1. The low-temperature susceptibility. 

Consider the low temperature series (IT^ for the full susceptibility [50] , for which 2043 
coefficients have been obtained in the u = variable j^D]. We denote this series 
F{u), so that 

F{u) = 41/2 -p 16 u® -k 104 -k 416 u® -k 2224 u® -k 8896 -k 43840 u® 

-k 175296 u® -k 825104 m®® -k • • • -k 02043 • + ••■ (24) 

Now consider this series modulo various integers q = 2”, {q = 2, 4, 8, 16, 32, 64, • • •) 
where we denote by Fq the corresponding series modulo q. We found the following 
simple results: 

F 2 {u) = O, Fi{u) = 0, F^{u) = Fieiu) = 4u^+8u\ (25) 

where the first two results are of no significance, and just reflect the lattice symmetry. 
However for <7 = 32 and q = 64, we found the appearance of simple lacunary series, 
so that 

F32{u) = 20 u2 -k24u^ + 16- ■ L{u), (26) 

Fg 4 {u) = 60u 2- (Il-k8u-kl0o2 + 8u^-k8u®) + (ASu'^ + 32u"^) ■ L{u), (27) 
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where L{u) corresponds to the first 1024 coefficients of the lacunary series with a 
natural boundary on the unit-circle |u| = 1 : 

n — oo 

■^(“) = 

n — 0 

= l+u + u'^ + U^ + u'^^ +U^'^ + -f -f -b -b • • • 

This strongly suggests that F 32 (u) and Fq^^u) satisfy the modulo 32 and modulo 64 
functional equations respectively: 

u^-F32iu) = F32iu^) -bl 6 u® -b24rt6 -bSrt®, (29) 

• (3 -b2u4) • F64 (u) = [3 +2u‘^) ■ Feiiu^) 

+ 16m®- (3 -b4M -b 6 u^ -b2M‘‘ -b58M® +4u^ + 2u'^ -b 6 u®®). (30) 

The series expansions Ti 28 , T 256 ) also satisfy similar functional equations, but 
they are more involved, the series having a less obvious lacunary series interpretation. 
For instance one finds that 

■ Fi 2 a{u) = Fi 2 a{u^) -b 32u® • (3 - u^) ■ L{u) -b 8 • u® • pis, 
where: P 13 = 2 — m — 8 — 15 m® -b 12 — 4 m® 

-b 8 M® -b4u^ - 8 u® -b4M“ - 8 u®®, (31) 

where L(v) is the lacunary series (|28)) . which satisfies the functional equation 
u + L{u^) = L{u). Therefore one deduces the functional equation modulo 128: 

u® • {u'^ - 3) • Fi28(m) - u'^ ■ {u^ - 2u^ - 3) ■ Fi2a{u^) 

= (m^ - 3) • Fi 28 ('u"‘) + 16 • P 28 where: (32) 

P28 = 4m 28_ 12^26 _ 2^24 _^gy 22 _^ 4^20 _ yl4_^ ^ 0^12 4^11 

— 2 -b 4 M® -b 12 U® — 10 — 13 M® — 11 U® -b 11 — 6 U® — — 3 u -b 3. 

Since we have seen that the full susceptibility series is quite close to the series 
expansion of x^‘^\ it is natural to ask if one obtains similar results modulo 2 ’", for 
From the series expansion (EU, we find that one obtains the same series as 
the one displayed in (1251) modulo 2, 4, 8 , 16. Modulo 32 and 64 one obtains simple 
functional equations for which are similar to (l 2 ^ and (IdHI) but actually slightly 
different. 

This can be rewritten in terms of the difference (1^^ . This difference (1^^ is zero 
modulo 2, 4, 8 , 16. Modulo 32 it is just one term, namely 16 u®® (the series for y^^^ 
being a non-trivial lacunary series) and modulo 64, it becomes the lacunary series 

Xl{v) - Xl\v) = 16m®® +32v'^° +32v^^ -b32M®® -b32M®® -b32M®®^ 

-b 32 M^®® -b 32 M®®® -b 32 m®° 2® -b 32 M^®®^ -b • • • (33) 

Remark: The low temperature series (IMl) relied on having coefficients up to 
the term in m^®'®®. Consequently the previous functional equations have been checked 
up to order 2043 in the the expansion (I24|) . The previous calculations underline the 
crucial role played by the lacunary series (l28)) where the next term is m^®®®. It would 
thus be interesting to validate a functional equation such as (13211 up to the point where 
the term m^®®® in (E51) is expected to emerge: this would require one to find jmsQ a few 


I Even though obtaining more terms for the low temperature series |[24t can be done with 
a polynomial time algorithm, getting more coefficients requires substantial computer resources: 
however, here the idea is that we just need a few extra terms. 
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(less than 10) extra terms in the low temperature series (1241) . Without trying to get 
more coefficients for the full susceptibility in exact arithmetic (not modulo a prime, 
or a power of a prime) which requires very substantial computer ressources, we can 
try to check all our previous functional equations modulo some integers of the form 
2’’, seizing the opportunity of having a polynomial algorithm to get many more than 
2000 coefficients (5000, 6000, 10000, ...) but just modulo 2, 4, 8, 16 ... 

3.2. The high-temperature susceptibility. 

Similarly, we now study the high-temperature expansion (1201) . modulo q with q = 
2, 4, 8, 16, 32, 64, and compare these series with the ones obtained modulo q with 
q = 2, 4, 8, 16, 32, 64 for (1^^ . Since, apart from the first constant coefficients, all 
the coefficients are divisible by 4, we introduce the series 

G(v) = ~ ^ = V +3v^ +8v^ +19v‘^ +Uv^ +I00v^ +224v'^ +491v^ 

4 

-f 1064 n® -f 2296 + 4932 + 10488 + 22180 -k • • • (34) 

and denote by Gq the corresponding series modulo q. We obtained the following 
results: Modulo 2 the series G 2 is the lacunary series L{v) — 1 (where L{v) is given 
by (I 2 HI)): 

G2{v) = v-i-v"^ -f -k n® -f -k -k -k -f -k -k • • • (35) 

which is a solution of the functional equation 

G2(v) = G2{v^) -\-v. (36) 

Modulo 4 the series G 4 is the lacunary series 3L{v) — 3 —2v 

G4(v) = V +3v^ -k3n^ +3v^ + 3 + 3 + 3 + 3 

+ 3v^^^ -k3u®^2 +3v^°^'^ + ■■■ (37) 

which is a solution of the functional equation: 

Giiv) = G^iv^) +v +2v^. (38) 

Modulo 8 the series Gs{v) becomes more difficult to recognise, 

Gsiv) = V + 30 "^ +3v^ + 3v^+ + + + lv^^ 

+ Av^'^ +Av^^ +Av'^° +Av^‘^ +Av‘^^ +Av^'^ +■■■ (39) 

though if we define 

Gs,{v) = Gs,{v) + L{v) -1, (40) 

then 

2- Gs,{v) = Av mod. 8. (41) 

Comparing the series (IMl) with the series (x^^ — l)/4 which is equal to 
1 / / 1 _ \ 

- • ( ( --^ ) — 1 ) = V -I- 3n^ -|- 8n® -I- 19 u"* -|- 44n® -I- 100 n® -I- 224 

4 VV(1 - 2z;)®/ / 

-f 490 n® -f 1064 w® -h 2296 + 4928 -k 10488 + 22180 -k • • • (42) 

one gets mod. 2, 4, 8, 16, 32 respectively: 

Vv'^v'^ mod. 2, z; -I-3 -I-3 -I-2 u® mod. A, 

z;-I-3 -I-3 -I-4 w®-I-4 u®-I-2 z;® mod. 8, (43) 

z;-k3n2-k8z;®-k3z;‘‘-kl2z;®-k4z;®-kl0z;®-k8z;^-k8ni°-k8z;i2-k8ni® mod. 16. 
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We see that, in contrast to the low-temperature expansion, a simple rational 
function like x^h yielding polynomial expressions modulo 2, 4, 8, 16, 32 cannot give 
rise to the emergence of lacunary series like (1351) and (I37p . For high-temperature series, 
one must therefore rather ask whether, modulo 2’’, one can distinguish between the 
full susceptibility xh and 

3.3. Functional equations mod. 2’’ for x the variable w. 

For the low and high-temperature series for x ia the variable w (see ([S]), ([7|). 
(0), (0, ...), we have obtained similar results and functional equations modulo 
q = 2, 4, 8, 16, 32, 64. These series, and the corresponding functional equations, 
are given in [Appendix A[ 

3.3.1. High temperature for x in the variable w. 

Let us consider the previous question of comparing xh with x*'^^ + but in the 
variable w, so that we are comparing xh and x^^^ + x*'^^ modulo 2^. 

The series expansion of the difference Ah = Xh — (x^^ + Xh'^) given in 
(O. The series expansion for x^^^ + X*'^^ can be obtained with an arbitrary number 
of exact coefficients, while 2043 coefficients of the series expansion of xh are known. 
Considering, modulo various integers, the 2043 coefficients of the series we found 
that Ah = 0 modulo 2’’, for r < 5. 

Modulo 16 one cannot distinguish xh and xiJ^ + X^^i their series expansions 
being a very simple lacunary series: 

Xh{w) = xiJ^ + ft^H — 10w+8w;^-|-8'u;® -I- Sic - L{w) (44) 

= 2w + +8w;® + 8w^'^ +8w^^ + 8w^^ + 8w^‘^^ + 8w‘^^'^ + 8w^^^ + ••• 

yielding the simple functional equation modulo 16: 

Xh{w'^) = w ■ XmIw) + 8w^ ■ {w^ — w —1). (45) 

Modulo 32, similarly, one cannot distinguish xh and x^^ + Xff\ their series 
expansions being a very simple lacunary series 

2w +8w^ +8w^ +24:w^'^ +24:w^^ +24w^^ + 24w^^^ + 24:w'^^'^ + 24:w^^^ + ■■■ 
= 10 ic -I- 16 -I- 8 + 8 ic® + 16 ic® + 24 lu • L{w), (46) 

where L{w) is the lacunary series (l28l) . This yields the simple functional equation 
modulo 32: 

Xni'uP') = w ■ xh{w) + 8w^- {2w^^ — + w —5). (47) 

Modulo 64, 128, similarly, one cannot distinguish between xh and x^^ + X^h + 
Xh'^, their series expansions being, again, very simple lacunary series. 

3.3.2. Low temperature for x in the variable w. 

/o\ 

Similarly, if one compares the low-temperature full susceptibility with xf (w) 
modulo 32 one finds the lacunary series: 

Xl(w) = 4ic^ -I-16ic® +24111® +16w;^^ +16i(;®® +16w;®® +16i(;®® 

+ 16 wi®® + 16 i(;®®° + 16 u;®i® + ••• , 


(48) 
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versus 

(w) = 4 + 16 + 24 w® + 16 + 16 + 16 + 16 w®® + 16 w®® 

+ 16 u>i®2 + 16 w^®® + 16 w®^® + 16 + ... ^ (49) 

the difference being only 16 w^®. 

One finds that the difference between xr and given in eqn. (I12L is zero 

modulo 2, 4, 8, 16, and equal to 16 ui^® modulo 32. Modulo 64 it is given by a 
lacunary series 

XL — Xl'^ = ■ L{w) +16w®® 


+ 32«;4. (tc^® -w^ -w - 


(50) 


One finds that the difference between xl and , as given by (TT^ , is zero 

modulo 2, 4, 8, 16, 32, 64 and is given by a lacunary series modulo 128: 

XL - {x^L + X^l) = 64'u;‘‘-L(w;) 

- +w^ +w^ +w^ +w +1). (51) 

If one includes yi®\ the difference (IT4l) between xl and + yi'^^ + y?^ is seen to 
be zero modulo 2, 4, 8, 16, 32, 64, 128, 256. 

The scenario seems to be that one cannot distinguish between the series for xl 
and that for a finite sum like y^^^ + y^^ + • • • + y^^"^ modulo 2’’ where r grows 
with n. The series expansion for the y^"^ are given in |23| . up to n = 12. This 
scenario has been checked and found to hold up to y)^ . Recall that the finite sum 
Xl ^ + Xl ^ + • • • + Xl is ^iso the diagonal of a rational function [TS] , implying that 
this finite sum reduces to algebraic functions modulo primes, or power of primes, and 
thus satisfies functional equations modulo primes, or power of primes. For instance, 
yL, which cannot be distinguished from this sum modulo 2’’ for some r, satisfies a 
functional equation modulo 2’’. 

These functional equations can thus be seen as related to the functional equations 

for y^^ For instance modulo 2, y^^(ru)/8 becomes the lacunary series 

~(3)/ \ 

Xh 


= w ■ Liw) — w ■ (ru^ + w‘‘' + ui +1), 

o 

from which one deduces the functional equation modulo 2: 

^. Xg M ^ xfjw'^) ^10 

8 8 

/o) 

Modulo 4 the series (w )/8 becomes the lacunary series 

~(3)/ N 

^ = 3w ■ L{w) + ic • (2 w® + +w‘^ + w +1), 

8 

from which one deduces the functional equation modulo 4: 

Xg (w) 


w ■ 
-14) 


X%\w‘^) 


+ W-{A+w^ - 2 w^®). 


(52) 


(53) 


(54) 


(55) 


For y))‘Vl6 we have similar results. The series Xl Vl® reduces, modulo 2, to 
■u;^®. Modulo 4 the series y^Vl® becomeslil the simple lacunary series 

~(4)/ \ 

Xl o4 7-/\| 16io4/28i 8, 4, 2, ii\ /r/?\ 

= 2w • L[w) -\-w ^-2w • [w w w +w -\-w-\-l), (56) 


16 

t Here, the calculations can be checked with an arbitrary number of coefficients. We did so with 
6000 coefficients. 
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yielding the functional equation modulo 4 


Xl\w^) 

16 


Xl M 

16 


+ • (2 


2 i (;16 + ^12 210 -^ - 1 ). 


(57) 


4. Automaton interpretation of the functional equations 

Recalling the decomposition of the full susceptibility as an infinite sum of n-fold 
integrals, these striking results can be seen as a consequence of the fact that, 

modulo integers that are powers of the prime 2, the full susceptibility series is the 
same lacunary series as the series for the first for instance the low-temperature 

series modulo 64 of the full susceptibility series and of -I- (which is the 
diagonal of a rational function) are the same. There is a not-widely-known discrete 
automaton [IIIIM] result that, modulo a prime p, diagonals of rational functions [38] 
not only reduce to algebraic functions, but also satisfy m “functional equations 
modulo p’'” of the form F{f{x), f{x^), ■■■ , f{x^ )) = 0. 

Let us recall some relevant results on discrete automaton |2Tl|24l|25l|26|: modulo 
a prime p, the diagonal of a rational function reduces to an algebraic function, and 
this is also true modulo p^ {p prime, r integer). Furthermore, these papers tell us 
that if f{x) is algebraic modulo a prime p then 1, f{x), P{x), p{x), - ■ ■ are linearly 
dependent, and 1, f{x), P{x), f^ (x), P {x)r" are also linearly dependent. From 
Fermat’s little theorem, namely that if p is a prime number, aP = a, (mod. p), one 
deduces for any series f{x) = ^ a„ ■ a;" 

(^a„-a;")P = '^al- x^^ = mod. p, (58) 

and, thus, f{xY = f{xP) modulo p, and, more generally, f{xY = /(x^ ) 

modulo p’'. One deduces that the relations F{f{x), /(x^), • • • , /(x^ )) = 0 can, in 
fact be written linearly, as 

X! PPP ■ /(^^") = 0, (59) 

n 

where the Pn(x) are polynomials with integer coefficients, (see for instance section 2 
in Lipshitz and van der Poorten [21jb 

Series generated by a finite automaton correspond to a system of algebraic 
eguations, which correspond, in turn (non trivially) to being algebraic. All these 
functional equations occurring for discrete automata can be seen as functional 
equations associated with algebraic functions modulo integers, in particular diagonals 
of rational functions. This can be seen as the origin of the functional equations of this 
paper. The functional equations we have obtained can be interprete43 as consequences 
of the fact that, modulo some integers that are powers of the prime 2, one cannot really 
make a distinction between the full susceptibility and the diagonal of a rational function 
(like the sum of the first x^"^s), and consequently reduce to algebraic functions modulo 

:( Equivalently, our conjectured functional equations can be seen as conjectures on the fact that, 
for instance, the non-holonomic infinite sum Xi — (x^l^ reduces to 

zero modulo 64, and possibly, that each series x^”^^ for n > 3, reduces to zero modulo 64, which 

corresponds to the (experimental) remark of section |[2j|, that the x^^\w) (resp. have an 

overall factor 2 '^. 
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2’’. One could say that non-holonomic functions, like the full susceptibility of the Ising 
model, correspond to “almost diagonal functions”. 

The automaton interpretation of this section can be revisited from a binomial 
viewpoint. Recall that the coefficients of the series expansion of diagonals of rational 
functions necessarily reduce to nested sums of products of binomials [271128] . Binomial 
coefficients modulo prime powers have been considered by many great mathematicians 
of the nineteenth centurjlj, yielding a large set of elegant results. Among the various 
prime powers, the powers of 2 seem to play a selectecITtl role [33]. Combining these 
two set of results is another approach to the main problem addressed in this paper, 
namely the study of (infinite sums of) diagonals of rational functions modulo prime 
powers. 

For powers of the prime 2, the functional equations satisfied by the full 
susceptibility are quite simple ones, which are associated with the lacunary seriesU (1351) . 
Of course for powers of other primes (S’", S’", ...), the functional equations satisfied 
by the full susceptibility should, if they exist, be much more involved, certainly not 
reducing to simple lacunary series. For powers of other primes, the scenario that 
modulo some powers of primes, one cannot disentangle the full susceptibility from 
some finite sum of is no longer valid. For instance, if one considers the series 

expansion (fTSl) of the difference xh — + Xh^)^ coefficient of 

and are, respectively, the following products of prime^ 

2 ^ 2 . 59 . 1403746269427- 1965616530023- 269691689798741092891, 

2® - 29 - 26811049 - 99658008281797903856656009433736710068597. (60) 

Similarly, if one considers the series expansion (|T3]) of one finds that the 

coefficient of and are, respectively, the following products of primes: 

2^-3^- 263- 3291604173673- 340864762033- 3935416959987419344918432619, 

2® - 5 - 5581 - 1400518348065785091954773485695960563962083118761649. (61) 

Besides powers of 2, there is no prime cancelling all the coefficients of the difference 
(fT8|) or of Xh^ . 

The question whether modulo primes different from 2, or powers of primes 
different from 2”, the full susceptibility, that no longer reduces to the sum of the 
first x^"^’s, satisfies (probably involved) functional equations remains open. 

5. Comments and speculations. 

5.1. Towards a physical interpretation of the functional equations. 

We can view these exact functional equations modulo integers that are powers of 
the prime 2, as a finite discrete automata [34] result corresponding to the fact that, 

J For instance Cauchy, Cayley, Gauss, Hensel, Hermite, Kummer, Legendre, see 12911301 . The study 
of congruences of combinatorial numbers ED usually starts with their p-adic order: it was first 
studied by Kummer m- 

tfln 1899 Glaisher observed that the number of odd entries in any given row of Pascal’s triangle is 
a power of 2. 

t Note that, as far as reduction to algebraic functions modulo powers of the prime 2 is concerned, 
a remarkably simple quadratic algebraic function corresponding to the Catalan number generating 
function also reduces to this lacunary series ll28t . as can be seen in [Appendix B| 

% Using the command “ifactor” in Maple. 
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modulo such integers, one cannot disentangle the full susceptibility from the diagonal 
of a rational function. From a more speculative, but more physical perspective, one 
might hope that such functional equations are the “shadow” modulo primes or powers 
of primes, of (probably very involved) functional equation^, the x ^ Frobenius 
symmetry being, in fact, an infinite order transformation. Furthermore, such an 
infinite order discrete transformation might be seen as a symmetry of the model. 
Along this symmetry line, recall that, in the case of the square Ising model, any 
isogeny of the elliptic curve parametrizing the model |10j can be interpreted as an 
exact generator of the renormalization group |35] . 

We remark that xf ' (see ([5])), can be written so that the Landen modulus clearly 
appears. Consider the 2 F 1 hypergeometric function $(a;), and recall the Landen 
modulus kL- 



From ([8]) x^l^ reads: 

xf = ^ik,) = 

= 4u"* +16u® +120u® + 480 ^ 1 ° +280011200^1^+ 58800 ^ 1 ® +••• 


One would like to see the Landen transformation k ^ kL, which can be viewed as 
an exact generate^ of the renormalization group [35] , as a symmetry of , before 
seeing it as a symmetry of the full susceptibility. Remarkably, this is the case: 4’(fci) 
and $(fc) are very closely and very simply related! This remarkable relation can be 
written in many ways, using the various variables we have introduced [g uni Eg El] 
(s, k, w, v), but since our functional relations are mostly written as series in v, we 
will write this relation in u. In u the Landen transformation corresponds to 




4 V 


1 + 4u2 


or: 


1 + 4u2' 


(63) 


Let us introduce 


4 /( 0 :) = 


1 + X (i<I>(a:) 


(64) 


X dx 

One then has the remarkably simple (differential-functional equation) relation 
representing the Landen transformation as a symmetry of : 

4 V 


*(TTb) = 


(65) 


Pursuing this line of argument on functional equations with an infinite order 
transformation (hopefully with a physical symmetry interpretation like the Landen 
transformation representing a generator of the renormalization group [35)1. it is 
tempting to imagine the u —^ infinite order transformation in functional equations 
like dMj), (El), or (I32|), as a mod. 2^ reduction of an infinite order symmetry of the 
model. In such a scenario, since one cannot distinguish, modulo 2 or 4, between v 


t In characteristic zero, not “modulo primes or powers of primes”. 

t This highly selected infinite order transformation (isogeny of the elliptic curve parametrizing the 
model [101135p has k = 1 as a fixed point, fc = 0, oo being clearly special ESj. 
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and v/{l + a functional equation G 2 {v) = G 2 {v^) + f, like (1551) . could also 
be written as 


G2 


( ^ \ 
V1 + 4 / 


G 2 {v^) + V or: 



( 66 ) 


and one could expect that the functional equations we discover modulo 2’’, are also the 
restriction modulo 2’' of some (quite involved) functional equations where the infinite 
order transformations have some physical meaning. Keeping in mind the unit circle 
natural boundary of the full susceptibility of the Ising model, it is worth recalling that 
functional equations like G{v) = G(u^)+u, not modulo integers but in characteristic 
zero, are the simplest examples to actually show that a series has a (unit circle) natural 
boundary. 

Recalling the expression of X^^\w) given by ([5]), the previous functional equation 
((65|) reads: 




VI + 4 


1 1 + 4ri^ 

8 dv 


(67) 


( 2) 

This might suggest replacing Xl {w) by xl{w) so as to consider this differential- 
functional equation (Ell) for the full susceptibility given by (El) modulo 2, 4, 8. 
Unfortunatly, in contrast with the calculations performed in section EH, one finds 
that a differentia^ functional equation like (1571) is not satisfied by the full susceptibility 
modulo 2, 4, 8. This seems to suggest that, if a “master” functional equation with an 
infinite order transformation symmetry exists (in exact arithmetic, not modulo some 
integers) for the full susceptibility, it is certainly much more involved that any simple 
generalization of (|67l) . 

Of course, all these ideas are quite speculative. It is reasonable to imagine 
that there must be some (probably involved) representation of the renormalization 
group [3S] of the full susceptibility. In particular, for this integrable model which 
has an elliptic parametrization, one might expect a representation of the action of the 
Landen transformation on the full susceptibility. 


5.2. Non-holonomic functions that are ratios of diagonals of rational functions, and 
beyond. 

Almost everything remains to be done to understand and describe this class of 
“nice” non-holonomic functions reducing to algebraic functions modulo some powers 
of primes. It is worth recalling that, while the product of two holonomic functions 
is a holonomic function, the ratio of two holonomic functions is, in genera^, non- 
holonomid The class of functions that are expressible as a ratio of two holonomic 
functions, and, further, the ratio of diagonals of rational functions, is clearly a very 
interesting class of functions: they are such that their series can be recast into series 
with integer coefficients [181138] (the ratio of series with integer coefficients is up to an 
overall integer a series with integer coefficients), and that their series, modulo primes, 
or modulo powers of primes, reduce to algebraic function^ (the ratio of series reducing 

X Therefore different from the functional equations in section 03}. 

f Except when the holonomic function in the denominator is an algebraic function: in that case the 
ratio is also holonomic. 

% More generally, a rational or even algebraic function (with integer coefficients) of holonomic 
functions, is such that it reduces modulo primes, or modulo powers of primes, to an algebraic function. 
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to algebraic functions reduces to the ratio of algebraic functions, and thus reduces to 
algebraic functions). Keeping in mind the (natural boundary of the) susceptibility of 
the Ising model, it is worth recalling that the ratio of holonomic functions can also 
yield a natural boundary, as can be seen from the solutions of non-linear Chazy III 
equations [3^l40] . 

The solutions of a particular non-linear third order differential equations having 
the Painleve property, the Chazy III equations [39l |40], have (circular) natural 
boundaries, and this can be seen as a direct consequence of the fact that the solutions 
correspond to the ratio of two holonomic functions, as shown by Chazy in crystal clear 
papers. 

The Chazy III equation |391140) is a third-order non-linear differential equation 
with a movable singularity that has a natural boundary for its solutions |41] : 



It can be rewritten in terms of a Schwarzian derivative: 

/(4) ^ 2f'^-{f,x} = 2 3/"2 with: y = ^- (69) 

Similarly, it is important to recall that the ratio of two holonomic functions, which is in 
general a non-holonomic function, is the solution of a non-linear Schwarzian derivative 
ODE: 

+R{x)- y = 0, t{x) = —, {t{x),x} = 2R{x). (70) 

dx 2/2 

The Chazy III non-linear differential equation (1681) has the quasi-modular form 
Eisenstein series E 2 / 2 . It can also be written as a log-derivativ^, namely a ratio 
A'/A, where Ramanujan’s modular discriminant function |42l 143) A is actually a 
selected holonomic function: a modular form. 

It is worth recalling, with the example of the enumeration of three-dimensional 
convex polygons |44) . that we have already encountered, in enumerative combinatorics, 
the emergence of ratios of holonomic functions. The class of functions characterised by 
ratios of holonomic functions and ratios of diagonals of rational functions is certainly 
an over-simplified scenario for the susceptibilitity of the Ising model. It is however 
an interesting “toy class” for the susceptibilitity of the Ising model, the class of 
the algebraic functions of diagonals of rational functions being much too large to 
reasonably explore. 

6. Conclusion 

This paper underlines the central role of discrete finite automata, or diagonals of 
rational functions, in lattice statistical mechanics and enumerative combinatorics, 
in particular regarding the challenging problem of the full susceptibility of the two- 
dimensional Ising model m 

The natural emergence of diagonals of rational functions in an extremely large 
set of lattice statistical mechanics and enumerative combinatorics models, has been 
emphasised and explained in m- That paper explains why a large class of functions 
describing lattice models that can be expressed as n-fold integrals of an algebraicQ 

I One takes the derivative with respect to the nome q (see equation (6) in |41|h 
% Or even holonomic. 
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integrand |18] . which are, consequently, solutions of linear differential equations, 
and, thus, at first sight, transcendental functions, is in fact a remarkable class of 
transcendental holonomic functions, namely diagonals of rational functions [11]. 

The corresponding selected linear differential operators are not only Fuchsian, but 
also |5] globally nilpotenlQ, and, since these transcendental functions are diagonals of 
rational functions, they reduce to algebraic functions modulo any prime [18]. They 
even reduce to algebraic functions modulo any integral power of a prime number. We 
may call this class of transcendental holonomic functions, that quite naturally occur 
in so many problems of theoretical physics [H] , “almost algebraic functions”. 

As far as transcendental non-holonomic functions are concerned, the full 
susceptibility of the two-dimensional Ising model is “algorithmically integrable” (with 
an 0{N^) polynomial algorithm) and can be decomposed as an infinite sum of n-fold 
integrals, that have been shown to be diagonals of rational functions [HI- Such nice 
transcendental non-holonomic functions emerging in physics require further concepts 
and tools to characterize and analyze them. 

In this paper we have obtained exact functional equations for low and high 
temperature series of the full susceptibility modulo integers that are powers of the 
prime 2, the series being associated with simple lacunary series. Since these 
exact results come from remarkably long low- and high-temperature series |20] with 
more than 2000 coefficients, these exact functional equations are currently not yet 
proved but extremely plausible conjectures. Recalling the decomposition of the full 
susceptibility as an infinite sum of n-fold integrals these striking results can, 

in fact, be seen as a consequence of the fact that, modulo integers that are powers 
of the prime 2, the full susceptibility series are the same series as the series for the 
sum of the first for instance the low-temperature series modulo 16 of the full 

susceptibility series and of (which is the diagonal of a rational function) are the 
same. 

Modulo a prime p, diagonals of a rational function not only reduce to algebraic 
functions, but also satisfy equations of the form F{f{x), f{xP), • • • , f[x^ )) = 0. In 
other words, the functional equations we have obtained, can be interpreted as the fact 
that modulo some integers that are powers of the prime 2, one cannot really distinguish 
between the full susceptibility and the diagonal of a rational function (like, for instance, 
^(2) -|-y(4)^ y scenario seems to be that one cannot distinguish the series for xl 
and for a finite sum like -I- x^^ -f • • • -f Xl modulo 2’' where r grows with n. 
The series expansion for the xi”^ £^re given in [23] . up to n = 12. This scenario can 
be checked up to Xl^^- Recall that the finite sum Xl^ + Xl^ -f • • • -f Xl"^ is also 
a diagonal of a rational function m, therefore this finite sum reduces to algebraic 
functions modulo powers of primes, and, thus, satisfies functional equations modulo 
powers of primes. Therefore xl which cannot be distinguished from this sum modulo 
some 2” satisfies a functional equation modulo some 2”. 

The question whether the full susceptibility satisfies (probably involved) 
functional equations modulo primes different from 2, or powers of primes different 
from 2”, remains open (even if, given ()60l) and (EB, it may seem unlikely). 

Much remains to be done to understand, and describe, this class of “nice” non- 
holonomic functions. It is worth recalling that, while the product of two holonomic 

f Their critical exponents are rational numbers, their Wronskian are A'^-th roots of rational functions, 
etc. 
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functions is a holonomic function, the ratio of two holonomic functions is, in general, 
non-holonomic. The class of functions that are expressible as ratios of diagonals 
of rational functions, is clearly a very interesting and important class of functions: 
they are such that their series (i) can be recast into series with integer coefficients, 
and (ii) modulo primes, or modulo powers of primes, reduce to algebraic function^ 
Concerning the susceptibility of the Ising model, it is worth recalling that ratios of 
holonomic functions can also yielc 0 a natural boundary. The ratio of diagonals of 
rational functions is probably an overly-simple scenario for the susceptibilitity of the 
Ising model. However it is clearly important to start studying this class of functions, 
and further, to study algebraic expressions of diagonals of rational functions, per se, 
before introducing them as a well-suited and powerful framework in which to study 
models of lattice statistical mechanics or enumerative combinatorics [H] . 
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Appendix A. Appendix: full susceptibility expansions in w 

Appendix A.l. Low-temperature expansion in w 

Let us consider the low-temperature expansion for xl in the w = \s/{l -\- s^) 
variable and introduce the series F{w) = xl/4: 

F(w) = ^ = -f350u>® -f 97020^^2 + 1585584 

4 

-f 25765744 -f 417159856 -f • • • (A.l) 

Modulo 2 and 4 the series (jA.ip becomes simple polynomials: 

F 2 {w) = mod. 2, Fi{w) = -|-2ri;® mod. 4. (A.2) 

Modulo 8 , this series becomes the lacunary series 
Fs{w) = w'^ -h6w® + F F ■■■ 

which satisfies the functional equation modulo 8 : 

F^(w^) — ■ Fs{w) F 2w^°- (2 Fw"^ — w^) = 0 mod. 8 . (A.3) 

Gomparing these results with x;^ j A, the series expansion (|S]) divided by 4 

= -f20w® -f350w® -h 5880^1° -f 97020 -f 1585584 

-f 25765740 -h 417159600 -f 6737127540 +•••, (A.4) 

f More generally, algebraic expressions of diagonal of rational functions are such that they reduce 
modulo primes, or modulo power of primes, to algebraic functions. 

% The fact that solutions of a particular Painleve-like non-linear third order differential equations, 
the Chazy III equations |39II4Q| . have (circular) natural boundaries is a direct consequence of the fact 
that the solutions correspond to the ratio of two holonomic functions, as shown by Chazy. 
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one finds that this series (IA.4I) gives, modulo 2, 4, 8, the same series expansions as 
(|A.2I) and (IA.3|) . and consequently satisfies the same functional equation as (IA.3I1 . 
Again, similarly to the results displayed in section (13.IL in the variable v, one cannot 


make, modulo 
-( 2 ) 

Xl ■ 


2, 4, 8, a distinction, for low-temperature expansions, between xl and 


Appendix A.2. High-temperature expansion in w 

Let us consider 0 , the high-temperature expansion for xh in the w = ^ 5 /( 1 -|-s^) 
variable, and introduce the series F{w) = xhIF. 

Fiw) = = w 4 16 -I-64?w^ -I- 256 w® -I- 1024 -|-4096 

-k 16384 w® -h 65540 -h 262144 -h 1048720 -k • • • (A.5) 

This series modulo 2 and 4 reads: 

F 2 {w) = w mod. 2, Fiiw) = w mod. 4, (A. 6 ) 

This series modulo 8 reads: 

Fs{w) = w +‘iw‘^ +4w^ 

-k -h 4w®^® -k ••• mod. 8 , 

which satisfies the functional equation: 

Fs{w‘^) -I- 4 w® • — w-I-1) = w ■ Fs{w) mod. 8 . (A.7) 

Modulo 16 it reads: 

Fiq{w) = w +4:w‘^ -k4u>® + 12 w^'^ -fl 2 w®® + 12 w^^ 

+ 12w^^^ +12w'^^’^ +12w^^^ + ■■■ mod. 16, (A. 8 ) 

from which one deduces the functional relation: 

Fie{w^) -I- 4u>® • {2 — w + 1) = w- Fiq{w) mod. 16. (A.9) 

Comparing the series (IA.5P with the series x^^/2, namely the series ([9|) divided by 2 
'( 1 ) 

= -—— = w +4,uP +16 uA + 6 ^X 0 ^^ -|-256w^ -I- 1024 w® 

2 1 - 4 w 

-k4096r(;^ -k 16384 w® -k 65536 w® -f 262144 -f 1048576 -f • • • (A.IO) 
one gets respectively, mod. 2, 4, 8 , 16, 32: 

w mod. 2, w mod. 4, w -|-4rt;^ mod. 8 , 

w -l-4'u;^ -|-16r(;® mod. 16, w -l-4'u;^ -|-16r(;® mod. 32, (A.11) 

w+4:w‘^+16uA mod. 6 A, w -|-4'u;^ -|-16 w® -I-64 mod. 128. 


Appendix B. A very simple algebraic function example illustrating the 
emergence of a lacunary series 


Let us consider a very simple algebraic function, the Catalan number generating 
function |33) : 


C{x) 


1 - VI -4a; 

2 X 


(B.l) 
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It is the solution of the quadratic equation x ■ C{xY — C{x) +1=0. Modulo 2 the 
series x ■ C(x) reduces to L{x) — 1, where L(x) is the lacunary series : 

L{x) = I + X + x^ + x"^ + x^ + x^^ + x^"^ + 

+ + x^°^^ + • • • (B.2) 

This Catalan number generating function (IB.II) satisfies the functional equation 

X ■ C(x^) — C{x) +1 = 0 mod. 2. (B.3) 

which can be seen as the consequence of C{x^) = C{x)‘^ mod. 2, or as the 

consequence of the functional equation L(pcf) = L(x) + x. 

This generating function (jB.ll) yields many other lacunary series modulo 2'’ (see 
for instance [53]), for instance, modulo 8 the series expansion of 4 + 4/(l —x- C{x'^)) 
reduces to 4 • L{x) where L{x) is given by (IB.21) . This result, namely the emergence 
of lacunary series, can be seen as a simple example of the previous finite automaton 
results, or, equivalently, congruence on algebraic functions, in the simplest case where 
only square roots occur. 
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